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Measurements of neutral pion (tv°) production at mid-rapidity in -^/sjvjv = 200 GeV Au + Au 
collisions as a function of transverse momentum, pr, collision centrality, and angle with respect to 
reaction plane are presented. The data represent the final n° results from the PHENIX experiment 
for the first RHIC Au + Au run at design center-of-mass-energy. They include additional data 
obtained using the PHENIX Level-2 trigger with more than a factor of three increase in statistics 
over previously published results for pr > 6 GeV/c. We evaluate the suppression in the yield of 
high-pT 7r°'s relative to point-like scaling expectations using the nuclear modification factor Raa- 
We present the pr dependence of Raa for nine bins in collision centrality. We separately integrate 
Raa over larger pr bins to show more precisely the centrality dependence of the high-pr suppression. 
We then evaluate the dependence of the high-pr suppression on the emission angle A(j> of the pions 
with respect to event reaction plane for 7 bins in collision centrality. We show that the yields of 
high-pr vr°'s vary strongly with A<f>, consistent with prior measurements We show that this 

variation persists in the most peripheral bin accessible in this analysis. For the peripheral bins we 
observe no suppression for neutral pions produced aligned with the reaction plane while the yield of 
7r°'s produced perpendicular to the reaction plane is suppressed by a factor of ~ 2. We analyze the 
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combined centrality and A(j> dependence of the n 1 suppression in different pt bins using different 
possible descriptions of parton energy loss dependence on jet path-length averages to determine 
whether a single geometric picture can explain the observed suppression pattern. 

PACS numbers: 25.75.Dw 
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I. INTRODUCTION 

High transverse momentum particles resulting from 
hard scatterings between incident partons have become 
one of the most effective tools for probing the proper- 
ties of the medium created in ultra-relativistic heavy ion 
collisions at RHIC. Data from the four RHIC experi- 
ments have unequivocally established the phenomenon 
of high transverse momentum hadron suppression in 
Au + Au compared to (appropriately scaled) p + p col- 
lisions [1, |3, H, H, 0> H, Q, while the lack of similar 
suppression in d + Au collisions [1, Qil HH, provides 
strong evidence that the suppression is not due to mod- 
ification of parton distributions in the incident nuclei. 
This suppression has been observed for a large variety of 
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hadron species, at highest pt for ir° and most recently 
rj 13] supporting further the notion that energy loss oc- 
curs at the parton level. Conversely, direct photon mea- 
surements by the PHEN1X collaboration show that the 
yield of hard photons in Au + Au collisions is consistent 
with p + p expectations scaled b y t he number of inco- 
herent nucleon-nucleon collisions [14j and, thus, provide 
final confirmation that hard scattering processes occur at 
rates expected from point-like processes. This observa- 
tion makes definitive the conclusion that the suppression 
of high-pr hadron production in Au + Au collisions is a 
final-state effect. Measurements of azimuthal angle cor- 
relations between hadron pairs resulting from fragmen- 
tation of hard-scattered partons into jets have provided 
additional confirmation of final-state medium effects on 
these partons [151 ]. 

Predictions of high-py suppression were made before 
the start of RHIC operation [n| and confirmation of 
these predictions may be considered one of the key 
successes of the RHIC program so far. The suppres- 
sion of high-pT single hadrons was predicted to result 
from the energy loss of hard-scattered quarks and glu- 
ons in the hot and dense QCD medium created in ultra- 
relativistic heavy ion collisions (see [It], Gj| an d refer- 
ences therein). In the canonical models, medium- induced 
gluon bremsstrahlung is expected to dominate the energy 
loss process and calculations of the high-p T sup- 

pression factor incorporating this effect have been able 
to successfully describe the experimental measurements 
[20l . l2ll I22I ]; however, recent measurements of heavy 
quark suppression pose some questions to this canoni- 
cal view. Nonetheless, from comparisons of the energy 
loss calculations with the experimental data, estimates 
of the initial net color charge density which is usually 
expressed in terms of a gluon rapidity density, dN g /dy, 
have been obtained yielding dN g /dy « 1000 and, assum- 
ing thermalization, estimates of the initial energ y d ensity 
have produced values in excess of 10 GeV/fm 3 [23l. [24|. 

However, in spite of this success, there are still a num- 
ber of outstanding issues with the interpretation of the 
Au + Au high-pT single-hadron suppression. Since the 
properties of the medium created in heavy ion collisions 
are not a priori known, the energy-loss calculations nec- 
essarily use the observed suppression to infer initial par- 
ton densities, usually through an intermediate param- 
eter that appears in the energy loss calculations. Al- 
though the initial parton density obtained by such "to- 
mographic" studies has to be consistent with the final 
(measured) total particle multiplicity, it is fair to ac- 
knowledge that the pr dependence of the suppression 
(rather than its absolute magnitude) is a more discrim- 
inating observable to test the various energy loss mod- 
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els. For tt° spectra, the suppression in central Au + Au 
collisions at ^/snn = 200 GeV is found to be approx- 
imately constant with px over the range, 3 < px < 
10 GeV/c. While the different energy loss calculations 
can reproduce this p-r-independent suppression, the de- 
tailed explanation of the constancy is different in each 
model. The effects invoked to explain the px depen- 
dence of the observed Au + Au high-py suppression in- 
clude: finite-energy effects, absorption of energy from the 
medium, evolution from incoherent (Bethe-Heitler) to co- 
herent (Landau-Pomeranchuk-Migdal or LPM) radiation 
with increasing parton energy [251 ]. the pr-dependent 
mixture of quark and gluon contributions to the hard- 
scattered parton spectrum, the increasingly larger ex- 
ponent of the underlying (power-law) parton px spec- 
tra [22j], and shadowing/EMC effect (26j. While most 
calculations of the high-p^ suppression in Au + Au col- 
lisions account for shadowing/EMC modifications of the 
nuclear parton distributions and for the relative mixture 
of quarks and gluons in the hard-scattered parton spec- 
tra, finite-energy corrections, absorption of energy from 
the medium, and the description of the energy loss pro- 
cess itself differs from calculation to calculation. Clearly 
the central Au + Au single-particle spectra are not suf- 
ficient, by themselves, to validate or exclude any of the 
different energy loss models; we must use more "differ- 
ential" probes of medium-induced energy loss to better 
understand the phenomenon. 

A robust prediction of non-Abelian parton energy loss 
calculations is that the average energy loss as a function 
of the in- medium path length L shows a quadratic depen- 
dence cx L 2 [13] ■ Such a behavior predicted for a static 
QCD medium, turns into an effective cx L-dependence 
in an expanding QGP 28]. In principle, the centrality 
dependence of the high-py suppression [g, @, Q provides 
an effective test of energy-loss calculations because the 
length of the path of the partons in the medium will 
change between peripheral and central collisions. How- 
ever, the energy loss calculations also have to account for 
changes in the initial properties of the medium with cen- 
trality and the extra flexibility in the description of the 
initial conditions means that the measured centrality de- 
pendence of the high-py suppression also does not strin- 
gently constrain energy loss models [29]. However, the 
path length of the parton in the medium can also be con- 
trolled by selecting high-py hadrons in different bins of 
azimuthal angle difference from the event-by-event deter- 
mined reaction plane. Indeed, shortly after experimen- 
tal observations of azimuthal anisotropy were reported 
[H EB] 1 arguments were made that the high-py anisotropy 
in non-central collisions was due to the spatial asymme- 
try of the medium and the resulting A(f> dependence of 
parton path lengths [3(1 ED]]. However, recent analyses 
have argued that the large azimuthal anisotropics at high 
Pt cannot be accounted for by energy loss alone - at 
least when realistic nuclear geometry is used to describe 
the spatial asymmetry of the initial state [1^, [32|, |33| . 
Some of these analyses were based on a picture of the 



energy loss process in which quarks or gluons that have 
emitted radiation effectively disappear from the steeply 
falling high-py spectrum because they are overwhelmed 
by partons of lower energy that escape from the medium 
losing little or no energy. In this picture, the medium 
effectively attenuates the high-p^ quarks and gluons and 
the high-py spectrum is dominated by partons originat- 
ing near the surface - i.e. partons originating in the 
"corona" [1^, [H, UK]. Then, the azimuthal anisotropy 
could be largely determined by the shape of the surface 
[3^| . However, it has been separately argued that fluctu- 
ations in the number of emitted gluons may be large and 
such fluctuations may weaken the corona effect [281 ] . 

In this paper we present measurements of n° produc- 
tion in ^snn = 200 GeV Au + Au collisions from the 
PHENIX experiment at RHIC. These data, obtained dur- 
ing Run-2 operation of RHIC in 2002, include additional 
data obtained with the PHENIX Level-2 trigger which 
improved the total statistics by a factor of ~ 3 com- 
pared to the prior analysis in [f| . The analyses presented 
here have also benefitted from advanced electromagnetic 
calorimeter calibrations and from improved understand- 
ing of the systematic errors in the 7r° measurement in 
course of the direct photon analysis presented in [Til ], 
where the ir° decay photons provide the main source 
of background. With the improved statistics, the px 
reach of the data is extended to higher px, allowing us 
to test whether the suppression starts to diminish above 
10 GeV/c in px- In addition, we extend the measure- 
ment of the centrality dependence of the suppression up 
to 8 GeV/c. 

We present measurements of the dependence of the tt° 
yield as a function of the angle A(f> of the 7r° with respect 
to the event reaction plane. By measuring the high-p^ 
hadron suppression as a function of A</>, for a given cen- 
trality bin, we can keep the properties of the medium 
fixed and vary only the average geometry of the jet prop- 
agation in the medium. By comparing different centrality 
bins we can, in principle, test how the initial properties of 
the medium affect the induced energy loss. Traditionally, 
measurements of the A<p dependence of hadron yields 
have been analyzed in terms of the elliptic flow param- 
eter, v 2 , and we note that the data presented here were 
used to obtain measurements of ir° v 2 for comparison to 
inclusive photon v 2 [Hj]. However, in this publication we 
focus not on v 2l but explicitly on the suppression as a 
function of Acf>, expressed in terms of the A(/)-dependent 
nuclear modification factor Raa(A(/>). While the data 
presented this way contain, in principle, the same infor- 
mation as the combination of A0-averaged Raa and v 2l 
Raa(A4>) provides a useful alternative way to evaluate 
the dependence of high-py suppression on geometry be- 
cause it effectively combines Raa(pt) and v 2 into a single 
set of data. We analyze the combined A</> and central- 
ity dependence of the high-p^ suppression in the context 
of different path-length and density dependencies of the 
parton energy loss process to evaluate whether any geo- 
metric picture can simultaneously describe the centrality 
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and Atfi dependence of the observed high pr deficit. 



II. EXPERIMENTAL DETAILS 

The data presented in this paper were obtained dur- 
ing Run-2 operation of the PHENIX experiment [35| at 
the Relativistic Heavy Ion Collider facility at Brookhaven 
National Laboratory (36|. The primary detectors used 
to obtain the presented results were the PHENIX cen- 
tral arm spectrometers, particularly the electromag- 
netic calorimeters (3?]], and the two beam-beam coun- 
ters (BBC's) (3l. In addition, the PHENIX zero-degree 
calorimeters [39j were used for triggering and centrality 
determination. 

Two-photon decays of neutral pions were measured in 
the PHENIX electromagnetic calorimeter, located at a 
radial distance of ~5.1 m from the beam-line, which 
has a pseudo-rapidity acceptance of —0.35 < 77 < 0.35 
and covers 7r radians in azimuth. The electromagnetic 
calorimeter is divided into eight sectors, with each sec- 
tor covering the full pseudo-rapidity range and ir/8 in 
azimuth. The calorimeter consists of two distinct parts 
using different technologies. A lead-scintillator sandwich 
calorimeter (PbSc) with 5 cm x 5 cm towers covers 3/4 
(6 sectors) of the central arm acceptance. A lead-glass 
Cerenkov calorimeter (PbGl) with 4 cm x 4 cm towers 
covers the remaining 1/4 (2 sectors) of the central arm 
acceptance. The corresponding A77 x A</> acceptance of a 
single tower at 77 = is 0.011 2 and 0.0075 2 for the PbSc 
and PbGl calorimeters, respectively. 

The event reaction plane in Au + Au collisions was 
measured in the two BBC's. Each BBC consists of 
64 hexagonal, quartz Cerenkov radiators closely packed 
around the beam pipe, in an approximately azimuthally 
symmetric configuration. The beam-beam counters, lo- 
cated 144 cm in each direction from the nominal center 
of the interaction diamond, are used to count charged 
particles produced in the pseudo-rapidity range 3.0 < 
< 3.9. The distribution of particles over the in- 
dividual channels of the BBC's allows measurement of 
the azimuthal distribution, dN^/ d(j>, of charged parti- 
cles within this pseudo-rapidity acceptance. The BBC's 
also provide measurement of the collision vertex position 
along the interaction diamond with a resolution of 0.6 cm 

m. 

The data presented here were obtained using the 
PHENIX minimum-bias Level-1 trigger, based on the 
BBC's and the PHENIX zero-degree calorimeters, that 
selects 92.2lg Q% of the total Au + Au hadronic interac- 
tion cross-section of 6.9 b @. For a subset of the data, 
events selected by the Level-1 trigger were subjected to 
software Level-2 trigger filtering after full assembly of 
events in the PHENIX event builder [4(|. A software 
algorithm performed a crude reconstruction of electro- 
magnetic clusters by summing the pedestal-subtracted 
and gain-calibrated energies of "tiles" made of adjacent 
4x4 calorimeter towers groups. The tiles are allowed 



to overlap such that every possible such tile that can be 
constructed in each calorimeter is tested. One of the 
Level-2 triggers ("LVL2A") selected events in which at 
least one cluster (tile) had energy > 3.5 GeV. Another 
Level-2 trigger ("LVL2B") selected events in the 50-92% 
centrality range (50% most peripheral events) with at 
least one cluster having energy > 1.5 GeV. 

The measurements presented in this paper were ob- 
tained from 31.4 M minimum bias triggers and approx- 
imately 1.7 M Level-2 trigger selected events. Of the 
Level-2 triggered events, 743 K events were selected by 
the higher energy LVL2A trigger and the remainder were 
selected by the peripheral, lower-energy LVL2B trigger. 
Taking into account their rejection factors, the two trig- 
gers sampled the equivalent of 44.4 M (LVL2A) and 
28.7 M (LVL2B) minimum-bias triggers, respectively. 
The difference is due to different online trigger pre-scale 
factors. Thus, the combined event sample contains ap- 
proximately a factor of 2.5-3 (considering both triggers 
over all centralities) more 7r°'s above 6 GeV/c than pre- 
viously published Run-2 ir° measurements [6J. 



III. DATA ANALYSIS 

A. Event Selection and Centrality 

In the offline analysis, the timing difference measured 
between the two PHENIX BBC's is used to determine 
the position of the collision vertex along the beam axis 
and to select events with vertex position within 30 cm of 
the nominal center of the detector for subsequent analy- 
sis. The energies measured in the zero-degree calorime- 
ters and the charged particle multiplicity measured in the 
BBC's are used to determine the collision centrality [4l| . 
For the ir° spectrum measurements presented here the 
total measured centrality range (0-92.2%) is subdivided 
into 9 bins: 0-10, 10-20, 20-30, 30-40, 40-50, 60-70, 70- 
80, 80-92.2%. For the reaction plane dependent analysis, 
the most central and two most peripheral bins are ex- 
cluded, the peripheral due to their large uncertainty in 
the reaction plane resolution, and the 0-10% bin simply 
because of its smaller intrinsic eccentricity. Additionally, 
we present also combined 0-20%, 20-60% and 60-92% 
data sets for comparison with other PHENIX analyses 
of high pt hadron production that use such centralities. 



B. Reaction Plane Measurement 

PHENIX has previously published measurements of 
elliptic flow using an event-by-event measured reaction 
plane [H, EH , and the same technique is used for the 
analysis presented here. Each BBC detector consists of 
128 quartz radiators placed in hexagonal, roughly con- 
centric rings whose light is individually collected by Pho- 
tomultiplier tubes (PMT's). The calibrated charge from 
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Centrality 

FIG. 1: Resolution of the reaction plane determined in the 
BBC versus centrality. The actual resolution is two times the 
square root of the quantity shown in the y-axis. 



TABLE I: Relative systematic uncertainty of the reaction 
plane resolution for the centrality bins used for the reaction 
plane analysis, (see Fig. [1} 



Centrality 



Syst. error 



0-10% 
10-20% 
20-30% 
30-40% 
40-50% 
50-60% 
60-70% 



20.3% 
5.1% 
3.9% 
3.8% 
4.1% 
4.6% 

22.5% 



each radiator is converted into an estimate for the num- 
ber of charged particles within the acceptance of each 
detector, iVj, using the measured single particle peak cen- 
troid. 

For the reaction plane measurement the measured N% 
values are corrected such that the weight of the inner 
rings which have the fewest radiators covering the full 
azimuthal angle range is reduced. Then, in terms of the 
corrected Ni values, \ the angle of the reaction plane 

is obtained from the formula 



tan (2*) 



E l ^ adj sin(20 l )-(E l iVf dj sin(2^)> 



E l ^f dj cos (24h) - (J2i N* a3 cos (200) 



adj 



(1) 

where represents the azimuthal angle of the center of 
a given radiator i. The subtraction of the average cen- 
troid position in Eq. {T]) removes the bias in the reaction 
plane measurement resulting from non-zero angle of the 
colliding beams, non- uniformities in detector acceptance, 
and other similar effects. The average is taken over many 
events localized in time with the event in question. A fi- 
nal correction is applied to remove non-uniformities at 
the 20% level in the if? distribution. 

Because the above-described procedure can also be ap- 
plied individually to each BBC, we have a redundant 



measurement of the reaction plane in the north and in 
the south, and we exploit this to determine the resolution 
of the full reaction plane measurement using standard 
procedures 44| . The resolution of the reaction plane is 
directly reflected in the quantity (cos 2 (^i — ^2)) where 
\L>i and ^2 are the reaction plane angles measured in each 
of the two beam-beam counters individually and the av- 
erage is taken over events. Figure Q] and Table. U show 
the variation of this quantity. The needed correction fac- 
tors can be derived from this using Eq. (fTTjl in Sec. IIII Fl 
i?^4^(A</)) Measurement, where the reaction plane correc- 
tions are described in more detail. 

The systematic errors associated with the measure- 
ment of the reaction plane come dominantly from how 
well the resolution is known. The uncertainty on this 
quantity is also shown with Figure [T] for all but the most 
peripheral centralities. This error is determined by ob- 
serving comparison of the similarly calculated quantity 
(sin2( , I'i — ^2)) which should by definition be equal to 
zero. The value of (sin2(^E' 1 — ^2)) is found to be con- 
sistent with for all centralities. The mean size of its 
fluctuations around are compared to the size of the 
(cos 2 (^i — ^2)) to derive the systematic errors in the 
table. Since the value of (cos 2 (^i — ^2)) decreases dra- 
matically in the lower multiplicity peripheral events, the 
relative size of the error increases. The size of this rel- 
ative error is also cross checked by comparing it to the 
relative error on elliptic flow (^2) measurements which 
is directly comparable since, as discussed in Sec. IIII Fl 
.Raa(A</>) Measurement, the resolution correction for Vi 
is a plain multiplicative factor. For the cross check, the v<x 
error is derived by taking the difference of t>2 made with 
reaction planes from the BBC North and BBC South 
separately. 

Because of the large rapidity gap between the PHENIX 
BBC's and the PHENIX Central Arm (An > 2.7 - 4.0), 
the measurements made in the BBC's are assumed to 
have no correlations (except collision geometry) with pro- 
cesses detected in the central arm that would affect the 
results presented in Sec. IIII F L Raa(A6) Measurement. 
Specifically, PYTHIA studies 45] indicate that any large 
rapidity-gap production correlated with jets (and thus 
the hard 7r°s we study) detected in the central arm have a 
negligible effect on reaction plane determination even for 
the most peripheral events considered in this paper. Fur- 
ther, we average both the North and South BBC, which 
are separated by An > 6.0, making potential effects of 
this nature especially unlikely. 



C. Neutral Pion Detection 

The detection of neutral pions is one of the major 
sources of information on identified particle production 
at high pr at RHIC, and PHENIX has already published 
the results of a number of ir° measurements in different 
colliding systems 0, H, [l(| 0, 53] • Here we will describe 
the technique for obtaining 7r° yields as a function of 
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Pt and centrality, which is now well established within 
PHENIX. 

Neutral pions are detected via their n° — > 7 + 7 de- 
cay channel. Due to the relatively short mean lifetime of 
neutral pions of about 10 -16 s, typical of electromagnetic 
decays, the pions decay close to the interaction point 
(ct rs 25 nm). This makes the decay vertex well known 
and the pions can be reconstructed via an invariant mass 
analysis of photon pairs measured by the EMCal. 

In the EMCal, hits or clusters are reconstructed by 
finding contiguous calorimeter towers with pulse heights 
above the ADC pedestal value. In order to obtain a 
cleaner sample of electromagnetic hits shower shape cuts 
are applied to select candidate photons and time-of-flight 
cuts are applied to reject slow hadrons. For the PbSc we 
require measured cluster times to be £ c iust < L/c± 1.2 ns 
where L is the straight-line path from the collision ver- 
tex to the reconstructed cluster centroid. For the PbGl 
we require reconstructed clusters to have times, i c iust < 
L/c ± 2 ns; the difference is due to the intrinsic timing 
resolutions of the two calorimeter technologies. 

The energy of each EMCal cluster is corrected for an- 
gular dependence and non-linearity based on test beam 
results and simulation. The linearity corrections for both 
detector types are different with the PbGl showing a 
stronger dependence on the energy. The correction fac- 
tors for a photon with a detected energy of 1 GeV (10 
GeV) are 1 (0.95) for the PbSc and 1.05 (0.975) for the 
PbGl, respectively The PbGl calorimeter also shows a 
stronger variation of the measured photon energy with 
the angle of incidence on the detector surface, at 20° the 
measured energy is reduced by 5% compared to perpen- 
dicular incidence (0°), while in the PbSc the effect is only 
of the order of 2%. 

In a typical Au + Au central event the EMCal de- 
tects > nearly 300 clusters corresponding to an occu- 
pancy of ~ 10% and therefore a non-negligible probabil- 
ity of cluster overlaps. To minimize the effects of cluster 
overlaps in high multiplicity events, the energy of each 
cluster in the PbSc calorimeter is determined not only 
from the sum of all contiguous towers with deposited en- 
ergy above a given threshold (15 MeV was our default 
value) but also, alternatively, "extrapolating" the mea- 
sured "core energy" of the 4-5 central towers assuming a 
standard electromagnetic shower profile in an event with 
zero background. For this latter case, the ecore energy 
was computed from the experimentally measured center 
of gravity, central shower energy, and impact angle in the 
calorimeter using a parameterized shower profile function 
obtained from electromagnetic showers measured in the 
beam tests. Such an ecore energy represented an esti- 
mate of the true energy of a photon impinging on the 
PbSc unbiased by background contributions from other 
particles produced in the same event and depositing en- 
ergy in the neighborhood of a given cluster. The use of 
ecore instead of the total cluster energy for photon re- 
construction, helped to reduce considerably the effects of 
cluster overlaps in central Au + Au collisions. 



For a photon pair originating from a ir° decay the in- 
variant mass 

m 77 = y \p it + P l2 f = ^2E 1 -E 2 ■ (l-cos0i 2 )(2) 

is identical to the ir° rest mass. However, due to the fi- 
nite energy and position resolution in the detection of the 
photon pair, the actual reconstructed value is smeared 
around a mean value, which can deviate from the nom- 
inal value. The reconstructed peak position is also in- 
fluenced by the high multiplicity in a heavy ion collision, 
where overlapping clusters can shift the measured energy 
of each photon. 

With the invariant mass analysis, 7r°'s cannot be iden- 
tified uniquely since all possible photon-photon combina- 
tions have to be considered. This leads to a large combi- 
natorial background, which increases quadratically with 
the multiplicity. The tt° yield is instead determined on a 
statistical basis, with the background contribution estab- 
lished via a mixed event technique as described below. 

One possibility to reduce the combinatorial back- 
ground is to make use of the phase-space distribution 
of the photons in a 7r° decay. For the ir° — * 7 + 7 decay, 
the two photons have minimum opening angle 

m 

tan6»i 2 /2 = — , (3) 
P 

where m is the 7r° mass and p its momentum, with p ~ p^ 
in the PHENIX central spectrometer. The angular dis- 
tribution of the 7 pair in the ir° rest frame, da/d cos 8*, 
is constant, which leads to a flat distribution in the mea- 
sured energy asymmetry of the two photons from ir° de- 
cay: 



where (3 = p/E ~ 1 is the velocity of the ir°. On the other 
hand, high px combinatorial pairs are strongly peaked 
near a = 1 because of the steeply falling spectrum of 
single photon candidates. This is illustrated in Fig. [2j 
where the asymmetry distribution for photons from 7r°s 
in a simulation is compared to the measured asymmetry 
for photon candidate pairs in real Au + Au collisions. In 
two independent analyses, asymmetry cuts of a < 0.7 
and a < 0.8 were employed, other values were used as a 
cross-check and to verify the energy scale (see below) . 

Pairs of candidate photon clusters within the PbGl 
and the PbSc calorimeter which satisfy the asymmetry 
cut are considered candidate 7r°'s. Fig. [3] shows exam- 
ple invariant mass distributions for ir° candidates with 
3.0 < pt < 3.5 GeV/c in Au + Au collisions for two dif- 
ferent bins of collision centrality. The background under 
the clear ir° mass peak in these figures is due to combina- 
torial mixing of photons from two different decaying 7r°'s 
or from pairs containing one or two non-photon clusters 
that nonetheless pass the above-described cuts. 

Such a combinatorial background can be determined 
by a so called mixed event technique. It is a widely 
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FIG. 2: Asymmetry of photon pairs with 3GeV/c < pr < 
5 GeV/c within the acceptance of one PbGl sector, for simu- 
lated single 7r°s and measured within minimum bias events. 
An asymmetry cut used during the analysis is also shown. 
(Due to the limited acceptance of the detector, the distribu- 
tion of the energy asymmetry shows a slight decrease towards 
a = 1). 



used method to determine the combinatorial background 
of combined particle properties, e.g. the invariant mass 
of a photon pair. The basic idea is to compare the re- 
sult obtained by combining particles within one event to 
the result for particle combinations from different events, 
which are a priori not correlated. 

In the case of the ir° invariant mass, the mixed event 
distribution is determined by combining one photon can- 
didate from the current event with all photon candidates 
from previous events. The number of previous events 
used for the pair combinations determines the statistical 
error of the background and is limited basically by com- 
puting resources. In this analysis a buffer of events is 
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FIG. 3: Invariant mass distributions of pairs of electromag- 
netic clusters passing photon selection cuts for pair trans- 
verse momenta satisfying 3.0 < pr < 3.5 GeV/c. Top panels: 
m 77 distributions in Au + Au events compared to a normal- 
ized mixed-event sample representing the combinatoric back- 
ground. Bottom panels: The m 77 distributions after subtrac- 
tion of the combinatoric background. Left: 10-20% centrality 
bin, Right: 60-70% centrality bin. 



varied from ~ 3-10 previous events is used for the event 
mixing with the current event. 

In order to describe the combinatorial background cor- 
rectly it is essential that the events used for mixing have 
similar properties as the real event and that they are 
not biased toward a certain reaction, as e.g. events that 
are chosen because of a high-px photon trigger. For this 
reason only minimum bias events are considered for the 
determination of the background in both data sets and 
different event classes for collision vertex, centrality and 
reaction plane are employed. 

It is self-evident that for the photons used in the event 
mixing the same criteria are applied as for the pair com- 
binations from one event, such as PID cuts, cuts on bad 
modules, and the asymmetry cut. Other properties valid 
a priori for the real photon pairs, e.g. a minimum dis- 
tance that allows to distinguish them, have to be consid- 
ered in addition. In the analysis a minimum distance cut 
of a least 8 cm is required for each photon pair combina- 
tion, within one event and for mixed events, respectively. 

For a given px bin the mixed-event background is nor- 
malized to the same-event invariant mass distribution 
outside the range of the ir° peak by scaling the mixed- 
event background with a function /(mi nv ). This scaling 
function is determined by fitting the ratio of the same- 
event and mixed-event invariant mass distribution for 
Pt bins up to 3GeV/c with a linear function. This is 
needed because at low px correlations in the real-event 
background due to overlapping clusters cannot be repro- 
duced by the mixed-event technique. For the px bins 
above 3GeV/c, a constant is used. To cross-check the 
result, the linear and the constant scaling function are 
also determined over the complete invariant mass region, 
including the 7r° peak, which is taken into account by an 
additional Gaussian in the fit function (e.g. a linear plus 
a Gaussian function). 

The determination of the scaling function for large 
pair-px is limited by statistics in the real event sample 
and does not lead to stable results. Instead a constant 
scaling factor is used if the ratio of the invariant mass dis- 
tribution shows bins with zero entries in the fit region. 
The scaling factor is determined by integrating the real 
and the mixed invariant mass distributions in the range 
with the peak region excluded. 

The scaled mixed-event background is subtracted from 
the same-event distribution to produce a statistical mea- 
sure of the true ir° yield. The result of such a subtraction 
procedure is shown in the bottom plots of Fig. [3J The 
raw 7T° yield is obtained in each px bin by integrating the 
subtracted invariant mass distribution in a range around 
the peak mean (m„o) of ± 3 times the Gaussian width 
(cr w o ) of the 7T° peak. Values of the mean and er^o , can be 
seen in Figure |4j Varying the size of the integration win- 
dow results in slightly different results, which contributes 
to the overall systematic uncertainty of the measurement, 
discussed in the Systematic Errors section MID 41 of the 
next chapter. 

Residual differences between the mixed background 
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and the foreground are still apparent in some px bins 
especially below ~ 2 GeV/c. Cluster merging, cluster 
splitting (fluctuations in the 2-D <f> — z energy profile 
cause multiple local maxima which are incorrectly sep- 
arated into distinct clusters), anti- neutron annihilation 
and even second order residual physics correlations such 
as 3 and multi-body decays, flow, HBT, etc., can all cause 
such differences. These remaining differences are com- 
pensated by the shape of the scaling function. In addi- 
tion, as a systematic check, the shape of the remaining 
background after subtraction is also fit with various low 
order polynomial functions and potential contributions 
to the peak yield are considered in the determination of 
the total systematic error from the peak extraction pro- 
cedure. 

The values of the peak width and mean are extracted 
in one initial analysis of the invariant mass distribution 
in which a p^-dcpcndcnt paramctrization is determined 
for different centralities. The use of predefined values for 
the position and spread of the tt° peak has the advantage 
that even in pt regions where no fit to the subtracted 
invariant mass distribution is possible, the integration 

I 



1. Acceptance and Detector Efficiency 

The geometric acceptance of the EMCal for the ir° — > 
77 decay is evaluated using a Monte-Carlo program that 
generates ir°s in a rapidity interval Ay with the same 
vertex distribution and rapidity distribution as observed 
in real events and contains the complete geometry infor- 
mation of the EMCal. The ir° decay is calculated via 
JETSET routines that are part of the PYTH1A event 
generator (48|. For each n° it is verified that both decay 
photons hit the detector. The resulting px distribution 
of accepted 7r°s is divided by the transverse momentum 
distribution of the generated 7r°s and provides the geo- 
metrical acceptance of the PbSc and PbGl, respectively. 

The detection efficiency is determined using GEANT 
to simulate the complete response of the calorimeter to 
single 7T° decays. The data from each simulated tt° is, 
then, embedded into real Au + Au events by adding the 
EMCal tower information of the simulated ir° to the 
tower information of the real event and recalculating the 
EMCal clusters. The efficiency for detecting the embed- 
ded 7T° is then again determined by comparing the input 
Pt spectrum to the reconstructed pr spectrum of the 
embedded 7r°'s. Using this technique we determine "effi- 
ciency" corrections that account for the energy resolution 
and position resolution of the calorimeter, as well as for 
the losses due to overlapping clusters in a real event envi- 



region is well defined just by extrapolation from low px- 



D. tv Spectrum Measurement 

For the reaction-plane independent ir° spectrum mea- 
surement in a given centrality class cent, the aforemen- 
tioned analysis is applied in Apx = 0.5 GeV/c bins for 
Pt > 1 GeV/c. We cease attempting to extract it yields 
at high px when the number of pairs within the selected 
(background-subtracted) 7r° mass window falls below 4 
counts. We then correct the resulting raw 7r° spectrum 
for the geometric acceptance a& y {pT), the overall detec- 
tion efficiency £ C ent(PT), which accounts for the cluster 
cut efficiency, the ir° mass cut efficiency, for losses due 
to cluster overlaps in high multiplicity events, for cuts on 
bad modules and for the calorimeter energy and position 
resolution. In addition a correction for conversion losses 
(cconv) in the material of the PHEN1X central arms and 
for the branching ratio of the two photon decay (c 77 ) is 
applied: 



(5) 

I 

ronmcnt. In addition, the embedding allows for a precise 
determination of the effect of edge cuts and bad mod- 
ules. Though these effects can be in principle considered 
as acceptance corrections, they depend not only on the 
geometry but also on the energy deposition of an electro- 
magnetic shower in the different calorimeter towers. 

In the embedding procedure the effects of photon con- 
versions are also included, as the GEANT simulation 
considers the material budget in front of the EMCal and 
the information for decay-photon conversions is retained. 
The final conversion correction, which is factorized from 
the rest of the efficiency for book-keeping purposes, is 
evaluated by comparing the 7r° yield with and without 
including conversions in the simulation. The final conver- 
sion correction, constant with pt depends on the photon 
PID cuts and material in front of each individual sec- 
tor and ranges from 6-8% in PbGl and 9-10% in PbSc. 
Comparing this to the sheer probability of a 7r° having 
at least one photon which converts, 21% PbGl and 14% 
PbSc, we see that a large portion of these tt° are still 
reconstructable. 

For the embedding, the input ir° spectrum is weighted 
to match a functional form fit to the measured tt° spec- 
trum so that the correct folding of the ir° spectrum with 
the resolution is obtained. This procedure is iterated, 
with the fit of the pt dependence of the input weights 
adjusted as the estimate of the efficiency correction im- 
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proves, until the procedure converges within the nearly 
p^-independent statistical error of the embedded sample, 
approximately 3%. 
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FIG. 4: Comparison of the 7r° peak position in real events 
measured with the PbGl and the PbSc with the results ob- 
tained for embedded 7r°s. The p-r-dependence of the peak 
position and width is due to calorimeter energy resolution as 
discussed in the text. 

Fig. 0] compares the invariant-mass peak after back- 
ground subtraction in the real data and the invariant 
mass peak of the embedded 7r° for the two different de- 
tector types. The measured n° peak position is shifted 
from the nominal value of approximately 134.98 MeV due 
to the finite energy resolution of the detector in combi- 
nation with the steeply falling spectrum and due to the 
additional effect of overlapping clusters. As illustrated 
the effects are well reproduced by the embedded 7r°'s. 



2. Trigger Efficiency 

The efficiency of the Level- 2 trigger is separately eval- 
uated by processing recorded minimum-bias events with 
the Level-2 trigger and evaluating the efficiency for the 
trigger to select events containing a high px cluster. This 
analysis shows complete (100%) efficiency for the LVL2A 
trigger at momenta well-above the trigger threshold of 
3.5 GeV/c (95% above 1.5 GeV/c for LVL2B) for obtain- 
ing clusters that also pass all offline cluster cuts. This 
is demonstrated in Figure [5] (a). The "plateau" values 
are determined from fitting the region above the turn-on 
also shown. 

The related trigger efficiency of reconstructed 7r°'s is 
calculated from fast MC based on these measured single 
cluster efficiencies. The calculation is performed both by 
using a integrated Gaussian fit to the single cluster effi- 
ciency and by directly using the finely binned histogram 



TABLE II: Corrections in the PbGl and PbSc to the raw tt° 
yield in central collisions (0-10%) and with TOF and shower 
shape cut applied. The main part of the efficiency loss in 
PbGl is due to the effect of bad module and edge cuts which 
is approximately 40% at high pr for the PbGl and 20% for 
the PbSc, respectively 





PbGl 


PbSc 




Pt 


3.25 GeV/c 8 


5 GeV/c 3.25 GeV/c 8 


.5 GeV/c 


dAi/=0.9 


0.068 


0.080 0.216 


0.246 




0.351 


0.358 0.455 


0.515 


Cconv 


0.93 


0.93 0.90 


0.90 






0.98798 





and constant plateau fit. Both methods give consistent 
results. The result for the latter method is shown in Fig- 
ure [5] (b), solid curves. The calculation is cross-checked, 
as demonstrated by the data points in Figure [5] (b) which 
show the ratio of the yield from the two Level-2 trigger 
samples per equivalent number of minimum bias events to 
the same from the true minimum bias sample itself. We 
combine the yields obtained in the minimum bias event 
sample and the LVL2A (LVL2B) trigger sample above 
a cut-off of 6.5 GeV/c (3.5 GeV/c) where the trigger 
reaches efficiencies greater than 0.4 such that the cor- 
rection factor is not allowed to be large. A conservative 
error of 3% is assigned to the efficiency calculations, re- 
sulting in a total error of ~ 3-5%, based on the three 
studies: 1) comparisons of the data shown in Figure El 
2) comparisons of the two calculational methods, and 3) 
a study of the yields in the subsample of minimum bias 
events which also fired the triggers, similar to 1). 



3. Other Corrections 

The calculated corrections are applied to the raw 7r° 
yield as given by Eq. ([5|) . Table [ill shows the correction 
in central collisions for two different bins in transverse 
momentum and for the PbGl and PbSc, respectively. As 
discussed above the effect of the cut on bad modules is in- 
cluded in the efficiency correction, due to its dependence 
on the depth of the electromagnetic shower. 

Following the usual PHENIX procedure of modifying 
the quoted yield values for each finite sized pt bin such 
that the measurement corresponds to pr value at the bin 
center instead of the average px of the bin [4!| (thereby 
facilitating taking ratios of spectra from different collision 
systems) , a final correction is applied to the yield of each 
data point. Using a continuous function which is fit to 
the data points, values for the invariant yields at the 
centers of the chosen pt bins are scaled by the ratio of 
the fit value at the fit's average pr to the fit value at the 
bin center. This is an iterative procedure similar to the 
final efficiency correction described above, with a smaller 
convergence criteria of < 0.1% of the previous correction. 
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FIG. 5: (color online) (a) Measured efficiency of single clus- 
ters of LVL2A (blue) and LVL2B (red) triggers as discussed in 
the text. The black lines are constant value fits to the plateau 
efficiency, greater than 99.7 (95%) for LVL2A (LVL2B). (b) 
Efficiency for neutral pion detection of the triggers as a func- 
tion of 7T° pr, calculated (solid curves) based upon the effi- 
ciencies in (a) and, as a cross check (data points) , compared 
to ratio of per equivalent minimum-bias event yields in the 
full trigger sample with the same in the true minimum bias 
sample. Since the latter is the ratio of two separate data 
samples, independent statistical fluctuations, as well as ~ 8% 
systematic effects in the yield extraction (discussed section 
IIIID 4[) in either sample can cause this measured ratio to be 
greater than 100%. 



4- Systematic Errors 

Each correction of the raw yield following Eq. (|5|) is af- 
flicted with its own uncertainty, but already the determi- 
nation of the ir° raw yield itself is sensitive to the method 
of extraction. In particular it is sensitive to the choice 
of the fit function for the background scaling and the ex- 
traction window. In principle, both should be taken into 
account by the detector efficiency, but in the efficiency 
calculation no background subtraction is necessary. For 
this reason the systematic error of the peak extraction 
method is determined in two steps: first via the compar- 
ison of the raw yield obtained with two different fits for 
the background scaling, and second through the compar- 
ison of the fully corrected spectra for different sizes of 
the extraction window, for the real data as well as for 



TABLE III: Summary of the dominant sources of systematic 
errors on the 7r° yields extracted independently with the PbGl 
and PbSc electromagnetic calorimeters in central events for 
different pr- For comparison the statistical uncertainty is also 
shown. 



PbGl 



PbSc 



p T (GeV/c) 


3.25 


8.5 


3.25 


8.5 


Yield extraction 


8.7% 


6.% 


9.8% 


7.3% 


Efficiency 


11.4% 


11.4% 


11.4% 


11.4% 


Acceptance 


2.5% 


2.5% 


2.5% 


2.5% 


Conversions 


2.9% 


2.9% 


2.9% 


2.9% 


Level-2 data 








3% 


Energy scale 


13.8% 


14.1% 


10.5% 


11.2% 


Total syst. 


20.5% 


19.3% 


18.8% 


18.7% 


Statistical 


10.6% 


50% 


8.1% 


26.6% 



the efficiency calculation. 

The systematic error introduced by the efficiency cal- 
culation is estimated by comparing the fully corrected 
spectra for different PID criteria as well as for different 
additional smearing. The smearing (or energy resolution 
in the simulation) is changed in a way that a clear dis- 
agreement between the measured 7r n peak width and the 
peak width from the embedding is observed. 

Apart from the uncertainty of the efficiency, the main 
contribution to the systematic error is the determination 
of the absolute energy scale. Based on the comparison of 
the 7T° peak positions in the data to the expectation from 
simulation the energy scale can only be determined or 
confirmed with limited accuracy, A(E)/E = 1.6% in the 
PbSc and, because of the smaller acceptance, A(E)/E = 
2% in the PbGl. 

The additional contributions to the systematic error 
that have not been discussed in detail involve the un- 
certainty of the conversion correction (2.9%) and of the 
acceptance calculation (2.5%) both due to small uncer- 
tainties in detector material and alignment. Table IIIII 
provides a final overview of the various contributions to 
the total error of the 7T° measurement in the PbSc and 
the PbGl, respectively. 

The most important cross check of the final result is 
the comparison of the result for the two different detec- 
tor types PbGl and PbSc, which is shown for peripheral 
events in Fig. [6j A good agreement within the errors 
is seen and similar consistency is found in all centrali- 
ties. Since they represent essentially independent mea- 
surements, the two results are averaged and the total 
error of the combined result is reduced using a standard 
weighted least-squares method also described in [5(3] • An 
additional cross check of the final result based on iso-spin 
symmetry is provided by the measurement of charged pi- 
ons in the central arm Q, this is shown for minimum- 
bias collisions in Fig. [7] The neutral pion measurement 
smoothly extends the result for charged pions to larger 
transverse momenta. 
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FIG. 6: Comparison of fully corrected spectra for the PbGl 
and the PbSc for peripheral events. Similar consistency is 
observed for all centralities. The error bars represent the sta- 
tistical and systematic uncertainties. 
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p T (GeV/c) 

FIG. 7: Comparison of the combined ty° result to the mea- 
surement of charged pions within the PHENIX experiment in 
minimum bias events. The fit considers the averaged result 
of the 7r + and ir~ measurement [3] below pr = 3GeV/c and 
the 7T° data above. 



E. R AA ( 



Measurement 



Using the invariant yields obtained from the above- 
described analysis and the separately measured invariant 
cross-section for ir° production in p + p collisions [46j , we 
calculate the nuclear modification factor, Raa, according 



to 



Raa(pi 



(l/N%"X)d*NZ A /dpTdy ^ 
(T A a) x d 2 a^°/dp T dy 1 



(6) 



where {Taa) is the average Glauber nuclear overlap func- 
tion for the centrality bin under consideration 



AA) 



f T AA (h)db 



j(i_ e --;r ,T ^( b ))db 



(7) 



from which the corresponding average number of 



nucleon-nucleon collisions, (N co u) 
easily obtained. 



r inel 



(T 



A A 



can be 



F. Raa(^-4>) Measurement 

The measurement of the raw 7r° yield with respect to 
the event reaction plane, A<fi = 4>(tt°) — VP, proceeds as 
described in Sec. IIII Dl for the px spectrum except that 
we measure the yields as a simultaneous function of both 
Pt and A<f>. Because the beam-beam counters have 2-7T 
acceptance, PHENIX can measure the tt° yields with uni- 
form acceptance over < A(j> < 2ir even though the 
electromagnetic calorimeters have only l7r nominal az- 
imuthal acceptance. Since the measurement of ^ is am- 
biguous wrt a 180° rotation of the reaction plane, and 
since we expect the ir° yields to be symmetric wrt reflec- 
tion around A<fi — 0, we measure the ir° yields in 6 bins 
of \Acf>\ over the range < \Acf>\ < ir/2. For each pt bin 
we evaluate the ratio, 



R(A<t>i,Pr) 



AN(A^p T ) 
E- =1 AiV(A^,p T ) 



(8) 



where N(A(f>^px) is the measured number of 7T°'s in a 
given (A(/),pt) bin, A<^ representing one orientation of 
Acj). Because the PHENIX central arm acceptance is ef- 
fectively constant as a function of A<fi and we do not ex- 
pect any azimuthal-dependence of our ir° efficiency cor- 
rections, R(A(f> i ,pT) can be written as: 



i?(A0 l5 p T ) = RAA{A(f>i,pT)/RAA(pT) 



(9) 



Using the measured Raa{pt) values we can directly 
convert the R(A(f>,pT) to Raa(A((),pt) without having 
to apply acceptance and efficiency corrections to the 
reaction-plane dependent yields. These corrections are 
already included in the Raa(pt) values as described 
above. 

However, before applying this procedure we must first 
correct the R(A<fi,pT) values for the finite resolution of 
the reaction plane measurement. One goal of our mea- 
surement is to determine Raa{A4>,pt) without assuming 
any particular functional dependence on A(f>. However, 
for purposes of correcting for reaction plane resolution, 
we take advantage of the fact that the observed ir° yields 
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FIG. 8: (color online) Neutral pion invariant yields as a func- 
tion of pr measured in minimum bias and 9 centrality classes 
in Au + Au collisions at ^snn ~ 200 GeV. Spectra for differ- 
ent centralities are scaled for clarity. Errors are total errors, 
full systematic and statistical added in quadrature. 
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FIG. 9: Nuclear modification factor Raa for neutral pions as 
a function of pr for different centralities. The shaded error 
band around unity indicates the uncertainty in scaling factor 
Taa and an overall scale uncertainty in the p + p reference. 



and hence the nuclear modification vary with A<p to first 
order as 



i? raw (A</>, pr) ~R {1 + 2v2 raw cos (2A0)) . 



(10) 



Where the superscript raw denotes the values not cor- 
rected for the reaction plane resolution. This resolution 
reduces v 2 by the factor ^2(^2 (W\ - * 2 )) 0, which 
is given by the independent measurement of \& in the two 
BBC's shown previously in FigQ]}). For each pr bin in 
a given centrality class we fit the R(A(j>,pr) values to 
the functional form in Eq. (|10p and then correct each 
measured R(A<P,pt) value according to 



R co "(A^ PT ) = R^(A ( j > ,p T ) 



1 + v 2 cmi cos (2A0) \ 
1 + w 2 raw cos (2A0) / 
(11) 



with v 2 co " = u 2 raw /V 2 ( cos2 (*i - *2)). We estimate 
the systematic error in the reaction plane resolution cor- 
rection by propagating the centrality dependent uncer- 
tainties in (cos 2 (^! — <3/ 2 ))from Fig QJ>). Of course, 
the above-described correction only strictly applies if 
Raa(A(/)) is well-described by the functional form in 
Eq. (fT0)l . While we do observe some departure from this 
harmonic form in the data, the differences are typically 
below 5% so our correction will not introduce a large 



IV. RESULTS 

A. 7T° Transverse Momentum Spectra and Nuclear 
Modification Factors 

The ir° invariant yields obtained using the procedure 
described in Sec. lIIIDl are presented in Fig. [8] as a func- 
tion of pr for the nine chosen centrality bins. With the 
increased statistics included in this analysis, we have ex- 
tended the pr range of the previous PHENIX measure- 
ment by at least 2 GeV/c for all centrality bins. The pr 
range of the central bin has been extended from 10 GeV/c 
to 14 GeV/c. Where the spectra overlap, the results 
shown here are consistent with the previously published 
results within systematic errors. The errors shown on 
the points in Fig. [5] include statistical errors and point- 
to-point varying systematic errors. The appendix tabu- 
lates the 7T° spectra plotted in Fig. [8] (centralities: 0-10%, 
10-20%, 70-80%, 80-92% ) plus the combined spec- 
tra for centralities 0-20%, 20-60% and 60-92%, which 
are used for comparison to other neutral meson mea- 
surements The spectra in Fig. [8] depart from the 
exponential-like shape above 3 GeV/c which is consistent 
with the expectation that high-p^ hadron production is 
dominated by hard-scattering processes which produce a 
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power-law px spectrum for hadrons resulting from quark 
and gluon fragmentation. 
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FIG. 10: Log-log plot of central Au + Au and scaled p + p it 
Pt distributions; 

In previous publications, we have established the sup- 
pression of high-pT 7T° production in Au + Au collisions 
[H, 0] . This suppression cannot be easily seen given the 
large range of invariant yield covered by Fig. [8] 

To evaluate the suppression of high-p^ 7r°'s, we show in 
Fig. [9] the px dependence of the 7r° nuclear modification 
factor, Raa(pt): f° r the nine individual bins of collision 
centrality and for the full minimum-bias centrality range 
0-92.2%. We make use of the PHENIX Run3 p + p base- 
line 7T° data. (47[. The error bars on the data in Fig. [9] in- 



clude contributions from statistical errors in the Au + Au 
and p + p measurements and from the systematic errors 
that do not cancel between the measurements. The sepa- 
rate band shown in each panel indicates p^-independcnt 
errors on the Raa measurement resulting from uncertain- 
ties in estimating Taa(x, y) and systematic errors on the 
normalization of the Au + Au and p + p measurements 
that do not cancel. As in previously published papers, 
{e.g. 0, H, 0]) a factor of ~ 5 high px tt° suppression 
in the most central Au + Au collisions, Raa ~ 0.2, is 
observed, with the suppression approximately px inde- 
pendent for px > 5 GeV/c. The suppression at high px 
decreases in more peripheral collisions such that the two 
most peripheral bins have Raa values consistent with 
unity for px > 3 GeV/c. 



B. Suppression via Spectrum Shift 

The suppression of high px particles as shown above 
was determined by comparison of the semi-inclusive mea- 
sured yields as a function of centrality in Au + Au colli- 
sions at y/sjvTV = 200 GeV to the (Taa) scaled px spec- 
trum from p + p collisions [5lj . A direct comparison of 
the 0-10% centrality Au + Au spectrum to the scaled 
p + p spectrum is shown in Fig. [TO] as a log-log plot to 
emphasize the pure power law dependence of the data 
for px > 3 GeV/c. The suppression is commonly ex- 
pressed by taking Raa the ratio of the point-like scaled 
semi- inclusive yield to the reference distribution (Eq. \§§ . 

As illustrated in Fig. [TOJ instead of viewing the sup- 
pression in the nuclear modification factor as "vertical" 
reduction of the Au+Au yields, it can equally well be 
taken as a "horizontal" shift in the (Taa) scaled Au + Au 
spectrum, such that 



J 



(l/N e /X)d 2 N A A(pT)/dp T dy 



(T 



AA 



d 2 Qpp{p T = Pt + S(p T )) 
dp' T dy 



r 



x (1 + dS(p T )/dpx) 



(12) 



where the last term in parenthesis is the Jacobian, 
dpx/dpx- 

Furthermore, owing to the pure power law of the 
p + p reference spectrum, Ed 3 cr/dp 3 cx p^ n with n = 
-8.10 ± 0.05 above px ~ 4 GeV/c, the relative shift of 
the spectra — assumed to be the result of energy loss for 
the Au + Au spectrum — is easily related to the equiva- 
lent ratio, Raa(j>t)' 



Raa(pi 



= (px + S( PT ))- n+1 



Pt 



n+l 



-(l + dS(px)/dp T ) (13) 



(1 + S( PT )/px)- n+i (l + dS( PT )/dpx) 



where the exponent is n — 1 because the relevant shift is 
in the da/dpx spectrum rather than in da/pxdpx- The 
fact that the Au + Au and reference p + p px spectra are 
parallel in Fig. [TO] provides a graphical illustration that 
the fractional px shift in the spectrum, S(px)/px = So, 
is a constant for all px > 3 GeV/c, which also results in a 
constant ratio of the spectra, Raa{pt)- For the constant 
fractional shift, the Jacobean is simply dS(px)/dpx = So 
and Eq. (fT"3|) becomes: 



Raa(pt) = (1 + S ) 



-n+2 



(14) 
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Raa(pt) 



l/(n-2) 



(15) 



The effective fractional energy loss, S\ oss , is related to 
the fractional shift in the measured spectrum, Sq. The 
hadrons that would have been produced in the reference 
p + p spectrum at transverse momentum pt + S(px) = 
(1 + Sq)pt, were detected with transverse momentum, 
Pt, implying a fractional energy loss: 



'los; 



1 - 1/(1 + S Q ) = 1 - R AA {p T ) 1,(n - 2) 



(16) 



The fractional energy loss S\ oss as a function of centrality 
expressed as -ZVp art is shown in Fig. [TT] for two different 
Pt ranges, 3 < pt < 5 GeV/c and 5 < pt < 7 GeV/c. 
There appears to be a small decrease of S\ oss with in- 
creasing pt, but the main observation from Fig.[TT]is that 

2/3 

S^ss increases approximately like N plirt , as suggested by 
GLV 52] and PQM [Ejjj]. 
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FIG. 11: (color online) Fractional energy loss S\ osa obtained 
from Eq. (|16p versus centrality given by N pa rt- The lines are 
fits of the form oc N^t f° r each pr range. 



It is important to realize that the effective fractional 
energy loss, S\ oss estimated from the shift in the pt spec- 
trum is actually less than the real average energy loss at a 
given pt, i.e. the observed particles havepr closer to the 
original value, than to the average. The effect is similar 
to that of "trigger bias" [54| where, due to the steeply 
falling spectrum, the (z) of detected single inclusive par- 
ticles is much larger than the (z) of jet fragmentation, 
where z — p^a ■ pj et /p 2 - et . Similarly for a given observed 

Pt, the events at larger p T with larger energy loss are 
lost under the events with smaller p T with smaller en- 
ergy loss. 

It should be noted that fluctuations due to the vari- 
ation of the path length and densities traversed by dif- 
ferent partons also contribute to the difference between 
<Sioss true and Sioss° bs - However, as long as the depen- 
dences of the induced energy loss on path length and 



parton energy approximately factorize, these fluctuations 
will also produce a py-independent reduction in Si oss ohs 
compared to S\ oss tlnc . 



C. Angle Dependence of High pr Suppression 

In order to try to separate the effects of the density 
of the medium and path length traversed, we study the 
dependence of the ir° yield with respect to the reaction 
plane. For a given centrality, variation of A</> gives a vari- 
ation of the path-length traversed for fixed initial condi- 
tions, while varying the centrality allows to determine 
the effect of varying the initial conditions. 
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FIG. 12: (color online) Raa vs. A(f> for 7r° yields integrated 
over 3 < pr < 5 GeV/c. Most statistical errors are smaller 
than the size of the points. The lines following the data points 
show the bin-to-bin errors resulting from the uncertainty in 
the reaction plane resolution correction (Fig. [1} and from bin- 
to-bin uncertainties in the Raa values. The shaded band 
indicates the overall Raa uncertainty. 

Figs. H2] and [T3] show the nuclear modification fac- 
tor Raa as a function of Ac/) integrated over pt £ 
]3GeV/c, 5GeV/c] and p T G ]5GeV/c, 8GeV/c], respec- 
tively. For all centralities (eccentricities) considered, 
there is almost a factor of two more suppression out- 
of-planc (Ac/) = tt/2) than in-plane (A<p — 0), something 
that is immediately apparent in viewing the data in this 
fashion, explicitly displaying information only implicit 
in Raa v 2, or the combination thereof. Strikingly, in 
contradiction to the data the variation in Raa with re- 
spect to the reaction plane expected by parton energy 
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FIG. 13: (color online) Raa vs. A(j> for tv° yields integrated 
over 5 < pr < 8 GeV/c. The error lines and band are the 
same as in Fig. 1121 



loss models [29|, [55[ should be much smaller for the more 
peripheral bins than for the central bins. As a result, 
the suppression vanishes (and perhaps an enhancement 
is observed) for smaller Acf> in the peripheral bins, corre- 
sponding to small path-length traversed in the medium. 
Although collective elliptic flow effects, usually not in- 
cluded in those models, are known to boost in-plane 
(compared to out-of-plane) particle production [2, |42| |. 
it is unclear how such collective effects can still play such 
an important role at the high pt bins considered. This 
may point to the possible need for a formation time be- 
fore suppression can occur [56j and which could also ex- 
plain why attempts to describe the azimuthal asymmetry 
V2 solely in terms of purely geometrical energy loss have 
failed. 

Figs. [14] and [15] give the angular dependence in terms 
of the fractional energy loss Si oss , and provide essentially 
the same information as shown in the plots of Raa{A-4>) 
in Figs. [T^l and [T31 Once again we see a large variation 
in energy loss as a function of angle. All the measure- 
ments of Raa or equivalently S\ oss vs. reaction plane 
and centrality, provide new constraints to models of jet 
quenching. To better understand the implications of the 
results shown in these figures, we will attempt in the next 
Section to find a common geometric description of the an- 
gle and centrality dependences in terms of an estimated 
path length of the parton in the medium. 



FIG. 14: (color online) Si oss vs. for tv° yields integrated 
over 3 < pt < 5 GeV/c. The statistical errors are smaller 
than the size of the points. The lines following the data points 
show the bin-to-bin errors resulting from the uncertainty in 
the reaction plane resolution correction (Fig. [1} and from bin- 
to-bin uncertainties in the S\ oss values. The shaded band 
indicates the overall S\ OS s uncertainty. 



D. Path-length Dependence of Energy Loss 

To analyze the path length dependence of parton en- 
ergy loss using the data presented here we will use dif- 
ferent methods for estimating the path lengths of par- 
tons in the medium as a function of centrality and 
Acf>. The "standard" approach would be to evaluate a 
length-weighted integral of the color-charge density in the 
medium along the parton path. We will adopt such an 
approach, described by the parameter L xy defined be- 
low, but we will also consider two other simplified ap- 
proaches that may help indicate which physics is most 
relevant in determining the observed suppression. We 
first consider, simply, L e , the distance from the edge to 
the center of the elliptical overlap zone of the Au + Au 
collision to represent the average path length of a parton 
in the medium. Then we try to weight the path length 
(or length-squared) traversed by a parton from the center 
of the ellipse by the participant density in the transverse 
plane: pL(A(f>) (pL 2 (Acf>)) . Finally, we do the same path 
length weighting for partons produced across the overlap 
ellipse, with hard-scattering production points weighted 
by TaaIx-.u)'- pL xy (pL 2 y ). It is obvious that such a 
A<p dependent analysis is not possible from just a simple 
combination of Raa and V2- 
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along the path length of the particle. This quan- 
tity is proportional to the opacity of the medium 
(n = L/X) divided by some undetermined cross- 
section. While the integral in principle extends to 
infinity the participant density naturally cuts off 
the integral outside the collision zone. 

/>oo 

pL= dr p paxt (r, A<j>). (17) 
Jo 

To account for the possible role of LPM coherence 
in the energy loss process, we evaluate a similar 
quantity including an extra factor of r in the inte- 
grand. 

/>oo 

P L 2 = / dr rp part (r,A0). (18) 
Jo 
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FIG. 15: (color online) S\ oas vs. Acj> for 7r° yields integrated 
over 5 < pr < 8 GeV/c. The error lines and band are the 
same as in Fig. 1141 

In detail, the three approaches considered here are as 
follows: 

(1) The simplest picture for the angular dependence of 

the energy loss in non-central collisions is that it 
is due to the asymmetric shape of the overlap re- 
gion of the colliding nuclei. Taking this idea to its 
extreme, only the simplest length scale, the length 
of the overlap region in a particular direction, mat- 
ters. To evaluate this length, we first estimated 
the RMS radius and eccentricity of an ellipse ap- 
proximating the shape of the overlap region from 
the transverse distribution of the participant den- 
sity calculated using standard Glauber techniques. 
We then estimated the path length, L e , of partons 
emitted at a given angle A</> by evaluating the dis- 
tance from the center of the ellipse to the edge. We 
thus ignore such considerations as event by event 
eccentricity fluctuations and only consider the av- 
erage. 

(2) Although the participant density is used to evaluate 

the dimensions of the ellipse, the above analysis 
ignores the dependence of participant density on 
position in the transverse plane. Thus as a natu- 
ral extension of the simple length scale in 1), for 
another analysis of the dependence of energy loss 
on A0, we assume that the color-charge density 
in the medium is proportional to participant den- 
sity and evaluate pL, the integral of this density 



We note that a Bjorken 1/r expansion of the 
medium would approximately cancel one power of 
r in the above expressions. Then, pL, might rep- 
resent LPM energy loss in the presence of 1-D ex- 
pansion. In the above integrals we assume all jets 
originate at the center of the collision region similar 
to our assumption for L e . 

(3) A final refinement on our geometrical calculation 
evaluates integrals like those in Eq. (fl~7| . (fT8| for 
jet production points distributed over the collision 
region to better account for geometric fluctuations. 
We are using a Monte-Carlo algorithm to sample jet 
production points (xo,yo) according to Taa{x,u) 
weighting and A</> angles from a uniform distribu- 
tion. For each jet, we evaluate the integral of the 
color-charge density (oc participant density) along 
the path of the parton out of the medium, 

pL xy = / dl p palt (xo + /cos A0, y + Ism A0). (19) 
Jo 

or 

/>oo 

pL 2 xy = / dl /ppart^o +/cosA0, y + / sin A0) . (20) 
Jo 

The above Monte-Carlo sampling yields a distribu- 
tion of pL xy (pL 2 xv ) values for each centrality. The 
larger values of pL xy (pL xy ) correspond to larger en- 
ergy loss which means these jets will have smaller 
contribution to the observed yield. To take this into 
account, a weighting factor is applied when evalu- 
ating (pL xy ). We assume that the energy loss can 
be represented by our empirical energy loss, S\ oss 
which we take to be proportional to pL xy (pL xy ) 
but with an undetermined multiplicative constant, 
K. We determine this constant in each centrality 
bin by relating S\ oss to Raa using Eq. p^|) and 
then evaluating the survival probability of each jet 
through 

P smv {pL xy ) = 1 - (k pL xy ) {n ~ 2) (21) 
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and requiring that the resulting suppression 
summed over all sampled jets agrees with the mea- 
sured A0-integrated Raa for that centrality bin. 
This determines the constant n(N part ) and allows 
us to evaluate a survival probability weighted aver- 
age for pL xy . 

We now evaluate how well the three above-described 
treatments of the geometry of the parton propagation in 
the medium perform in providing a consistent description 
of the A0 and centrality dependence of ir° suppression. 
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FIG. 16: (color online) Raa and Sioss vs. L e whose definition 
is explained in the text. Each data point represents a central- 
ity bin and A<f> (azimuth defined w.r.t. the reaction plane) bin 
combination. The six centrality bins are denoted by different 
colors as follows: cyan 60-70%, mauve 50-60%, blue 40-50%, 
green 30-40%, red 20-30%, black 0-10%. Within each central- 
ity group, the six different data points correspond to the same 
A(j> bins as in Figs. 1121 - 1151 The height of the bars around 
each data point represent the systematic error in .Raa(A^) 
(iSioaa corresponding to L e ). 

The plots shown in Figs. [1511201 illustrate the path 
length and color-charge dependence of suppression us- 
ing our empirical estimators. The systematic error in 
the estimators due to the uncertainty of the overlap ge- 
ometry parameter in a centrality class are approximately 
10-20% and is not shown in the figures. This uncertainty 
is derived by propagating the impact parameter and ec- 
centricit y u ncertainties from the PHENIX Glauber MC 
itself P,gl|. 

From Figs. [T51l2T)l it is evident that the individual cen- 
trality bins exhibit roughly parallel linear dependencies 
of the variables vs pL xy etc. For the 3 < pt < 5 GeV/c 
bin these slopes are such that the curves are disjoint due 
to the steeper value of the slopes in each centrality group 
(each color in the plots) compared to the bin-to-bin trend. 
For the higher 5 < px < 8 GeV/c bin, the slopes in the 
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FIG. 17: (color online) Raa and Si oas vs. pL, the density 
weighted path-length. Colors/data points as in Fig. 1161 
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FIG. 18: (color online) Raa and Si oss vs. pL 2 , the density 
weighted path-length squared. Colors/data points as in Fig. 



individual centralities flatten such that they follow the 
bin-to-bin trend much better. These are meant to be 
qualitative statements. We defer anymore quantitative 
tests, e.g. statistical tests, to subsequent datasets (e.g. 
the larger PHENIX 2004-2005 Run4 dataset) with which 
we can improve statistical precision. 

In this spirit, we note several other interesting quali- 
tative dependencies: 
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FIG. 19: (color online) Raa and S\ OS s vs. pL xy whose defi- 
nition is explained in the text. Colors/data points as in Fig. 

m 



L e <2 fm. 

• Neither Raa nor S\ oss is universal as a function of 
pL, pi? or pL xy for 3 < pt < 5 GeV/c. 

• For the higher 5 < px < 8 GeV/c px bin, Si OS s 
(itU^) approaches universality as a function of pL 2 , 
pL, and pL xy (possibly to a lesser extent for the lat- 
ter two) but does not achieve the level of universal- 
ity found for L £ . The largest deviations from uni- 
versality in these quantities are towards the longer 
axis (perpendicular to the event plane) in the more 
peripheral events. The dependence of S\ oss is rea- 
sonably linear as a function of pL but tends to level 
off at larger values of pL 2 . 

• When pL is normalized by the central density 
P P art(0,0) = p C ont, then S\ oss (Raa) become uni- 
versal in the quantities pL/ p ccnt for both px ranges 
with a linear dependence. The universality appears 
to become more exact in the higher px range. A 
similar improvement (not shown) of the qualitative 
universality for pL 2 and pL xy is also observed when 
these quantities are scaled in the same way by p C ent- 
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FIG. 20: (color online) Raa and Si 08S vs. pL normalized by 
the most central {b x = b y = 0) density p C cnt- Colors/data 
points as in Fig. [TS] 



• Raa is universal as a function of L s for all centrality 
classes and both px ranges. 

• S'loss is universal and is a linear function of L e for 
all centrality classes and both px ranges. 

• Within our errors, we see no suppression Raa ~ 1, 
hence no apparent fractional energy loss S\ ass for 



The most important of these observations is the ab- 
sence of suppression for the same value of L e < 2 fm for 
both px ranges, 3 < px < 5 GeV/c and 5 < px < 8 
GeV/c. This may suggest a "formation time effect" (see 
[Ifili [13) or some other type of emission zone which has 
generally not been taken into account in parton energy- 
loss models. The level of universal scaling with this sim- 
ple geometric quantity is surprising. 



V. SUMMARY AND CONCLUSIONS 

We have presented a detailed analysis of high-pr neu- 
tral pion suppression as a function of transverse momen- 
tum, centrality and angle with respect to reaction plane 
in Au + Au collisions at ^/snn = 200 GeV. The tt° yields 
have been measured in the range px ~ 1 - 14 GeV/c in 
nine centrality bins and compared to the ir° differential 
cross-sections measured in p + p . The ratio of Au + Au 
over p + p spectra (scaled by the number of equivalent 
nucleon-nucleon scatterings) is reduced more and more 
for larger centralities. The resulting suppression factor 
is, however, independent of px above px ~ 4 GeV/c for 
all centralities. This observation can be interpreted as an 
indication of a constant effective fractional energy loss, 
fixed S'loss "pr shift" , in the Au + Au compared to the 
p + p yields. The dependence of S'loss in the centrality as 
given by the number of participating nucleons N part fol- 

2/3 

lows an approximately N palt law as predicted by parton 
energy loss models. 

In order to constrain the "jet quenching" models with 
more differential observables, we have experimentally 
tested the path-length (L) dependence of the energy loss 
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by exploiting the spatial azimuthal asymmetry of the sys- 
tem produced in non-central nuclear collisions. Due to 
the characteristic almond-like shape of the overlapping 
matter produced in A+A reactions with finite impact pa- 
rameter, partons traversing the produced medium along 
the direction perpendicular to the reaction plane ("out- 
of-plane") will comparatively go through more matter 
than those going parallel to it ( "in- plane" ) , and there- 
fore are expected to lose more energy. 

We have studied the suppression pattern along differ- 
ent A(j) trajectories with respect to the reaction plane 
determined with the Beam-Beam-Counters at high ra- 
pidities. The measured Raa(A<P) curves show clearly a 
factor of ~2 more suppression out-of-plane (A(p = tt/2) 
than in-plane (A<fi — 0) for all the centralities (eccen- 
tricities) considered. Theoretical calculations of parton 
energy loss in an azimuthally asymmetric medium pre- 
dict a significantly smaller difference between the sup- 
pression patterns for partons emitted at A</> = and 
Acj) = tt/2 [2^, [H, UK]. The discrepancy is stronger for 
more peripheral centralities (with correspondingly larger 
eccentricities) and challenges the underlying in-medium 
path-length dependence of non-Abelian parton energy 
loss. Although elliptic flow effects are responsible for 
extra boost of in-plane (compared to out-of-plane) pions, 
it is unclear how such collective effects persist up to pr 
values as high as ~ 8 GeV/c. We have analyzed the ob- 
served reaction-plane and centrality dependence of the 
nuclear modification factor with three different versions 
of a Monte Carlo model with an increasing level of refine- 
ment in the description of the azimuthal propagation of 
the parton in the medium. For all three approaches we 
observe that the ir° suppression tends to vanish for values 
of the path-length L sa 2 fm in the two pt ranges consid- 
ered, 3 < p T < 5 GeV/c and 5 < p T < 8 GeV/c. Such a 
result suggests either a formation time effect or a surface 
emission zone that results in a pT-hidependent suppres- 
sion and puts additional constraints to parton energy-loss 
models. 
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spectra 



TABLE IV: Final combined PbSc+PbGl tt invariant yields 
vs. pr for centrality 0-10%. 



PT 


Yield 


Stat. Error 


% 


Sys. Error 


% 


1.25 


3.314E+00 


2.518E-02 


0.76 


4.026E-01 


12.15 


1.75 


5.981E-01 


4.946E-03 


0.83 


6.784E-02 


11.34 


2.25 


1.208E-01 


1.253E-03 


1.04 


1.447E-02 


11.98 


2.75 


2.718E-02 


3.744E-04 


1.38 


3.521E-03 


12.96 


3.25 


6.970E-03 


1.270E-04 


1.82 


9.751E-04 


13.99 


3.75 


2.158E-03 


4.713E-05 


2.18 


2.686E-04 


12.44 


4.25 


7.185E-04 


2.133E-05 


2.97 


9.349E-05 


13.01 


4.75 


2.715E-04 


1.063E-05 


3.92 


3.575E-05 


13.17 


5.25 


1.288E-04 


5.931E-06 


4.61 


1.702E-05 


13.21 


5.75 


5.417E-05 


2.606E-06 


4.81 


7.731E-06 


14.27 


6.25 


2.940E-05 


1.560E-06 


5.31 


4.106E-06 


13.97 


6.75 


1.280E-05 


9.501E-07 


7.43 


1.922E-06 


15.02 


7.25 


7.641E-06 


6.459E-07 


8.45 


1.241E-06 


16.24 


7.75 


4.630E-06 


4.668E-07 


10.08 


7.508E-07 


16.22 


8.50 


1.883E-06 


1.809E-07 


9.61 


3.033E-07 


16.11 


9.50 


1.057E-06 


1.276E-07 


12.07 


1.952E-07 


18.47 


11.00 


2.777E-07 


4.274E-08 


15.39 


5.664E-08 


20.39 


13.00 


5.941E-08 


1.704E-08 


28.87 


1.222E-08 


20.57 
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TABLE V: Final combined PbSc+PbGl ir° invariant yields 
vs. p T for centrality 10-20%. 



Pt 


Yield 


Stat. Error 


% 


Sys. Error 


% 


1.25 


2.054E+00 


1.461E-02 


0.71 


2.655E-01 


12.93 


1.75 


4.137E-01 


2.933E-03 


0.71 


4.616E-02 


11.16 


2.25 


8.576E-02 


7.654E-04 


0.89 


1.039E-02 


12.11 


2.75 


2.028E-02 


2.305E-04 


1.14 


2.612E-03 


12.88 


3.25 


5.057E-03 


7.980E-05 


1.58 


6.778E-04 


13.40 


3.75 


1.665E-03 


3.170E-05 


1.90 


1.995E-04 


11.98 


4.25 


5.859E-04 


1.511E-05 


2.58 


7.301E-05 


12.46 


4.75 


2.253E-04 


7.948E-06 


3.53 


3.003E-05 


13.33 


5.25 


9.486E-05 


4.369E-06 


4.61 


1.246E-05 


13.14 


5.75 


4.651E-05 


2.087E-06 


4.49 


6.696E-06 


14.40 


6.25 


2.224E-05 


1.249E-06 


5.62 


3.252E-06 


14.62 


6.75 


1.109E-05 


8.621E-07 


7.78 


1.899E-06 


17.13 


7.25 


6.455E-06 


5.485E-07 


8.50 


1.091E-06 


16.90 


7.75 


3.568E-06 


3.999E-07 


11.21 


7.173E-07 


20.10 


8.50 


1.724E-06 


1.718E-07 


9.96 


3.279E-07 


19.01 


9.50 


6.318E-07 


9.789E-08 


15.49 


1.144E-07 


18.11 


11.00 


1.701E-07 


3.347E-08 


19.68 


3.147E-08 


18.51 


13.00 


5.093E-08 


1.610E-08 


31.62 


9.747E-09 


19.14 



TABLE VIII: Final combined PbSc+PbGl tt° invariant yields 
vs. pr for centrality 40-50%. 



Pt 


Yield 




Stat. Error 


07 
70 


— 

Sys. Error 


07 

70 


1.25 


6.389E-01 


3.367E-03 


0.53 


7.216E-02 


11.29 


1.75 


1.156E-01 


6.789E-04 


0.59 


1.315E-02 


11.37 


2.25 


2.442E-02 


1.926E-04 


0.79 


2.911E-03 


11.92 


2.75 


6.172E-03 


6.521E-05 


1.06 


7.890E-04 


12.78 


3.25 


1.682E-03 


2.455E-05 


1.46 


2.194E-04 


13.04 


3.75 


5.822E-04 


1.161E-05 


1.99 


7.179E-05 


12.33 


4.25 


1.927E-04 


6.113E-06 


3.17 


2.480E-05 


12.87 


4.75 


8.818E-05 


o a T r* T7i f\r- 

3.476E-06 


3.94 


1.169E-05 


13.26 


5.25 


3.627E-05 


2.166E-06 


5.97 


4.995E-06 


13.77 


5.75 


1.611E-05 


9.656E-07 


5.99 


2.261E-06 


14.04 


6.25 


9.635E-06 


6.880E-07 


7.14 


1.490E-06 


15.47 


6.75 


4.467E-06 


4.278E-07 


9.58 


7.232E-07 


16.19 


7.25 


2.044E-06 


2.585E-07 


12.65 


3.197E-07 


15.64 


7.75 


1.363E-06 


2.198E-07 


16.13 


2.882E-07 


21.15 


8.50 


7.878E-07 


1.056E-07 


13.41 


1.409E-07 


17.88 


9.50 


2.197E-07 


5.630E-08 


25.62 


4.969E-08 


22.61 


11.00 


1.053E-07 


2.280E-08 


21.66 


2.116E-08 


20.10 


13.00 


2.792E-08 


1.140E-08 


40.82 


6.121E-09 


21.92 



TABLE VI: Final combined PbSc+PbGl tt° invariant yields 
vs. pt for centrality 20-30%. 





Pt 


Yield 


Stat. Error 


% 


Sys. Error 


% 



1.25 


1.601E+00 


9.668E-03 


0.60 


1.852E-01 


11.57 


1.75 


2.879E-01 


1.911E-03 


0.66 


3.260E-02 


11.32 


2.25 


6.045E-02 


5.117E-04 


0.85 


7.416E-03 


12.27 


2.75 


1.429E-02 


1.537E-04 


1.08 


1.761E-03 


12.32 


3.25 


3.983E-03 


5.534E-05 


1.39 


5.192E-04 


13.04 


3.75 


1.233E-03 


2.340E-05 


1.90 


1.546E-04 


12.53 


4.25 


4.749E-04 


1.158E-05 


2.44 


6.115E-05 


12.88 


4.75 


1.732E-04 


5.898E-06 


3.41 


2.258E-05 


13.04 


5.25 


7.761E-05 


3.503E-06 


4.51 


1.074E-05 


13.84 


5.75 


3.573E-05 


1.627E-06 


4.55 


4.870E-06 


13.63 


6.25 


1.714E-05 


9.568E-07 


5.58 


2.389E-06 


13.94 


6.75 


9.015E-06 


6.625E-07 


7.35 


1.384E-06 


15.36 


7.25 


5.146E-06 


4.423E-07 


8.59 


8.214E-07 


15.96 


7.75 


2.878E-06 


3.267E-07 


11.35 


5.465E-07 


18.99 


8.50 


1.363E-06 


1.452E-07 


10.65 


2.517E-07 


18.46 


9.50 


6.216E-07 


8.347E-08 


13.43 


1.088E-07 


17.50 


11.00 


1.825E-07 


2.972E-08 


16.28 


3.299E-08 


18.08 


13.00 


3.552E-08 


1.267E-08 


35.68 


6.852E-09 


19.29 



TABLE VII: Final combined PbSc+PbGl tt° invariant yields 
vs. pr for centrality 30-40%. 





Pt 


Yield 


Stat. Error 


% 


Sys. Error 


% 



1.25 


1.040E+00 


5 


.648E-03 


0.54 


1 


244E-01 


11. 


96 


1.75 


1.754E-01 


1 


100E-03 


0.63 


2 


001E-02 


11 


41 


2.25 


3.833E-02 


3 


102E-04 


0.81 


4 


567E-03 


11. 


91 


2.75 


9.610E-03 


9 


930E-05 


1.03 


1 


.175E-03 


12. 


23 


3.25 


2.670E-03 


3 


.764E-05 


1.41 


3 


.512E-04 


13. 


15 


3.75 


8.612E-04 


1 


.667E-05 


1.94 


1 


097E-04 


12. 


74 


4.25 


3.270E-04 


8 


.158E-06 


2.49 


4 


.185E-05 


12. 


80 


4.75 


1.252E-04 


4 


421E-06 


3.53 


1 


619E-05 


12. 


94 


5.25 


5.266E-05 


2 


.822E-06 


5.36 


7 


394E-06 


14. 


04 


5.75 


2.761E-05 


1 


.348E-06 


4.88 


3 


.839E-06 


13 


90 


6.25 


1.189E-05 


8 


.138E-07 


6.85 


1 


.949E-06 


16. 


39 


6.75 


7.115E-06 


5 


804E-07 


8.16 


1 


.198E-06 


16. 


84 


7.25 


3.705E-06 


3 


.972E-07 


10.72 


6 


264E-07 


16. 


91 


7.75 


1.898E-06 


2 


.549E-07 


13.42 


3 


307E-07 


17. 


42 


8.50 


1.168E-06 


1 


301E-07 


11.13 


1 


967E-07 


16. 


83 


9.50 


5.043E-07 


8 


.312E-08 


16.48 


9 


.634E-08 


19 


10 


11.00 


1.541E-07 


2 


.748E-08 


17.83 


2 


910E-08 


18. 


89 


13.00 


2.941E-08 


1 


.278E-08 


33.46 


5 


.621E-09 


19 
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TABLE IX: Final combined PbSc+PbGl tt° invariant yields 
vs. pt for centrality 50-60%. For points with no errors given, 
data value represents 90% confidence level upper limit. 



Pt 


Yield 


Stat. Error 


% 


Sys. Error 


% 


1.25 


3.593E-01 


1.941E-03 


0.54 


4.022E-02 


11.19 


1.75 


6.197E-02 


4.018E-04 


0.65 


7.069E-03 


11.41 


2.25 


1.309E-02 


1.175E-04 


0.90 


1.553E-03 


11.87 


2.75 


3.479E-03 


4.211E-05 


1.21 


4.205E-04 


12.09 


3.25 


1.019E-03 


1.695E-05 


1.66 


1.291E-04 


12.67 


3.75 


3.480E-04 


8.518E-06 


2.45 


4.380E-05 


12.59 


4.25 


1.329E-04 


4.558E-06 


3.43 


1.763E-05 


13.26 


4.75 


4.959E-05 


2.434E-06 


4.91 


6.310E-06 


12.73 


5.25 


2.125E-05 


1.585E-06 


7.46 


3.032E-06 


14.27 


5.75 


9.917E-06 


7.569E-07 


7.63 


1.540E-06 


15.52 


6.25 


6.127E-06 


5.471E-07 


8.93 


9.978E-07 


16.29 


6.75 


3.246E-06 


3.392E-07 


10.45 


4.965E-07 


15.30 


7.25 


1.664E-06 


2.449E-07 


14.72 


3.102E-07 


18.65 


7.75 


1.129E-06 


1.886E-07 


16.70 


2.114E-07 


18.72 


8.50 


3.362E-07 


7.419E-08 


22.07 


6.694E-08 


19.91 


9.50 


1.817E-07 


4.619E-08 


25.42 


3.329E-08 


18.32 


11.00 


2.858E-08 


1.112E-08 


38.89 


4.803E-09 


16.81 


13.00 


2.311E-08 












TABLE X: Final combined PbSc+PbGl 


7T° invariant yields 


vs. pt 


for centrality 60-70%. 










PT 


Yield 


Stat. Error 


% 


Sys. Error 


% 



1.25 


1.731E-01 


1.121E-03 


0.65 


1.985E-02 


11.47 


1.75 


3.022E-02 


2.288E-04 


0.76 


3.425E-03 


11.33 


2.25 


6.567E-03 


7.011E-05 


1.07 


7.773E-04 


11.84 


2.75 


1.644E-03 


2.565E-05 


1.56 


2.057E-04 


12.51 


3.25 


5.255E-04 


1.158E-05 


2.20 


6.682E-05 


12.72 


3.75 


1.801E-04 


6.044E-06 


3.36 


2.259E-05 


12.54 


4.25 


6.986E-05 


3.184E-06 


4.56 


9.254E-06 


13.25 


4.75 


2.312E-05 


1.631E-06 


7.06 


3.101E-06 


13.41 


5.25 


1.156E-05 


1.145E-06 


9.90 


1.720E-06 


14.87 


5.75 


4.884E-06 


5.045E-07 


10.33 


7.560E-07 


15.48 


6.25 


2.690E-06 


3.650E-07 


13.57 


4.303E-07 


16.00 


6.75 


1.822E-06 


2.658E-07 


14.58 


3.369E-07 


18.48 


7.25 


6.281E-07 


1.480E-07 


23.57 


1.178E-07 


18.76 


7.75 


2.446E-07 


1.082E-07 


44.22 


4.632E-08 


18.94 


8.50 


1.417E-07 


4.482E-08 


31.62 


2.707E-08 


19.10 


9.50 


1.094E-07 


3.843E-08 


35.14 


2.106E-08 


19.26 


11.00 


2.492E-08 


1.114E-08 


44.72 


4.816E-09 


19.33 


13.00 


4.728E-09 


4.728E-09 


100.00 


9.226E-10 


19.51 
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TABLE XI: Final combined PbSc+PbGl tt° invariant yields TABLE XIV: tt° spectrum for combined centralities: 0-20% 



vs. pt 


for centrality 


70-80%. 








Pt 


Yield 


Stat. Error 


% 


Sys. Error 


% 


Pt 


Yield 


Stat. Error 


% 


Sys. Error 


% 


1.25 


2.684E+00 


1.455E-02 


0.54 


3.106E-01 


11.58 


1.25 


7.416E-02 


5.166E-04 


0.70 


8.842E-03 


11.92 


1.75 


5.059E-01 


2.875E-03 


0.57 


5.131E-02 


10.14 


1.75 


1.282E-02 


1.189E-04 


0.93 


1.496E-03 


11.67 


2.25 


1.033E-01 


7.343E-04 


0.71 


1.137E-02 


11.01 


2.25 


2.721E-03 


3.774E-05 


1.39 


3.245E-04 


11.92 


2.75 


2.373E-02 


2.198E-04 


0.93 


2.837E-03 


11.96 


2.75 


7.455E-04 


1.514E-05 


2.03 


9.131E-05 


12.25 


3.25 


6.014E-03 


7.501E-05 


1.25 


7.693E-04 


12.79 


3.25 


2.461E-04 


7.508E-06 


3.05 


3.248E-05 


13.20 


3.75 


1.912E-03 


2.840E-05 


1.49 


2.137E-04 


11.18 


3.75 


7.200E-05 


3.689E-06 


5.12 


9.687E-06 


13.46 


4.25 


6.522E-04 


1.307E-05 


2.00 


7.660E-05 


11.74 


4.25 


2.609E-05 


2.071E-06 


7.94 


4.034E-06 


15.46 


4.75 


2.484E-04 


6.637E-06 


2.67 


3.055E-05 


12.30 


4.75 


1.288E-05 


1.308E-06 


10.15 


2.161E-06 


16.78 


5.25 


1.118E-04 


3.683E-06 


3.29 


1.368E-05 


12.23 


5.25 


4.650E-06 


7.727E-07 


16.62 


9.050E-07 


19.46 


5.75 


5.034E-05 


1.670E-06 


3.32 


6.775E-06 


13.46 


5.75 


2.416E-06 


3.897E-07 


16.13 


4.736E-07 


19.60 


6.25 


2.582E-05 


9.994E-07 


3.87 


3.466E-06 


13.43 


6.25 


1.763E-06 


2.713E-07 


15.39 


2.795E-07 


15.85 


6.75 


1.194E-05 


6.415E-07 


5.37 


1.825E-06 


15.29 


6.75 


5.945E-07 


1.651E-07 


27.77 


1.221E-07 


20.53 


7.25 


7.048E-06 


4.237E-07 


6.01 


1.115E-06 


15.82 


7.25 


4.817E-07 


1.245E-07 


25.84 


8.088E-08 


16.79 


7.75 


4.099E-06 


3.073E-07 


7.50 


7.159E-07 


17.46 


7.75 


1.344E-07 


6.718E-08 


50.00 


2.545E-08 


18.94 


8.50 


1.804E-06 


1.247E-07 


6.92 


3.037E-07 


16.84 


8.50 


1.135E-07 


4.012E-08 


35.36 


2.167E-08 


19.10 


9.50 


8.445E-07 


8.042E-08 


9.52 


1.491E-07 


17.65 


9.50 


4.968E-08 


2.484E-08 


50.00 


9.568E-09 


19.26 


11.00 


2.239E-07 


2.714E-08 


12.12 


4.219E-08 


18.84 


11.00 


5.060E-09 


5.060E-09 


100.00 


9.778E-10 


19.33 


13.00 


5.517E-08 


1.011E-08 


18.32 


1.044E-08 


18.93 



TABLE XII: Final combined PbSc+PbGl ^° invariant yields TABLE XV: tt spectrum for combined centralities: 20-60% 
vs. pr for centrality 80-92%. 















Pt 


Yield 


Stat. Error 


% 


Sys. Error 


% 


PT 


Yield 


Stat. Error 


% 


Sys. Error 


% 


1.25 


9.097E-01 


2.963E-03 


0.33 


9.041E-02 


9.94 


1.25 


3.494E-02 


6.093E-04 


1.74 


4.504E-03 


12.89 


1.75 


1.602E-01 


5.854E-04 


0.37 


1.570E-02 


9.80 


1.75 


6.037E-03 


1.291E-04 


2.14 


7.607E-04 


12.60 


2.25 


3.407E-02 


1.599E-04 


0.47 


3.579E-03 


10.51 


2.25 


1.319E-03 


3.628E-05 


2.75 


1.701E-04 


12.89 


2.75 


8.386E-03 


4.969E-05 


0.59 


9.138E-04 


10.90 


2.75 


3.321E-04 


1.243E-05 


3.74 


4.570E-05 


13.76 


3.25 


2.339E-03 


1.832E-05 


0.78 


2.711E-04 


11.59 


3.25 


1.059E-04 


5.281E-06 


4.99 


1.483E-05 


14.01 


3.75 


7.562E-04 


8.034E-06 


1.06 


8.387E-05 


11.09 


3.75 


3.625E-05 


2.408E-06 


6.64 


4.455E-06 


12.29 


4.25 


2.819E-04 


4.021E-06 


1.43 


3.259E-05 


11.56 


4.25 


1.233E-05 


1.293E-06 


10.48 


1.730E-06 


14.03 


4.75 


1.090E-04 


2.126E-06 


1.95 


1.263E-05 


11.58 


4.75 


6.501E-06 


7.988E-07 


12.29 


9.044E-07 


13.91 


5.25 


4.695E-05 


1.310E-06 


2.79 


5.959E-06 


12.69 


5.25 


3.018E-06 


5.360E-07 


17.76 


4.224E-07 


13.99 


5.75 


2.234E-05 


6.107E-07 


2.73 


2.907E-06 


13.01 


5.75 


1.072E-06 


2.315E-07 


21.60 


1.815E-07 


16.94 


6.25 


1.120E-05 


3.833E-07 


3.42 


1.605E-06 


14.33 


6.25 


3.265E-07 


1.154E-07 


35.36 


5.945E-08 


18.21 


6.75 


5.961E-06 


2.591E-07 


4.35 


8.809E-07 


14.78 


6.75 


2.805E-07 


9.918E-08 


35.36 


5.185E-08 


18.48 


7.25 


3.140E-06 


1.732E-07 


5.52 


4.937E-07 


15.73 


7.25 


2.231E-07 


8.434E-08 


37.80 


4.187E-08 


18.76 


7.75 


1.817E-06 


1.264E-07 


6.96 


3.292E-07 


18.12 


7.75 


8.467E-08 


4.888E-08 


57.74 


1.604E-08 


18.94 


8.50 


9.139E-07 


5.846E-08 


6.40 


1.581E-07 


17.30 


8.50 


3.602E-08 


2.080E-08 


57.74 


6.880E-09 


19.10 


9.50 


3.818E-07 


3.462E-08 


9.07 


7.055E-08 


18.48 


9.50 


1.077E-08 


1.077E-08 


100.00 


2.074E-09 


19.26 


11.00 


1.176E-07 


1.194E-08 


10.15 


2.061E-08 


17.52 


11.00 


4.375E-09 


4.375E-09 


100.00 


8.455E-10 


19.32 


13.00 


2.899E-08 


6.031E-09 


20.80 


4.497E-09 


15.51 



TABLE XIII: Final combined PbSc+PbGl tt° invariant yields 
vs. pr for centrality 0-92%. 





PT 


Yield 


Stat. Error 


% 


Sys. Error 


% 



1.25 


1.078E+00 


3.333E-03 


0.31 


1.205E-01 


11.17 


1.75 


1.928E-01 


6.847E-04 


0.36 


2.171E-02 


11.26 


2.25 


4.038E-02 


1.742E-04 


0.43 


4.822E-03 


11.94 


2.75 


9.578E-03 


5.293E-05 


0.55 


1.202E-03 


12.55 


3.25 


2.564E-03 


1.858E-05 


0.72 


3.375E-04 


13.17 


3.75 


8.115E-04 


7.353E-06 


0.91 


1.013E-04 


12.48 


4.25 


2.906E-04 


3.475E-06 


1.20 


3.729E-05 


12.84 


4.75 


1.121E-04 


1.806E-06 


1.61 


1.466E-05 


13.08 


5.25 


4.924E-05 


1.031E-06 


2.09 


6.494E-06 


13.19 


5.75 


2.240E-05 


4.723E-07 


2.11 


3.012E-06 


13.45 


6.25 


1.190E-05 


2.909E-07 


2.44 


1.647E-06 


13.83 


6.75 


5.970E-06 


1.943E-07 


3.25 


8.494E-07 


14.23 


7.25 


3.246E-06 


1.273E-07 


3.92 


4.758E-07 


14.65 


7.75 


1.715E-06 


9.049E-08 


5.28 


2.658E-07 


15.49 


8.50 


8.583E-07 


3.892E-08 


4.53 


1.285E-07 


14.98 


9.50 


3.078E-07 


2.351E-08 


7.64 


5.041E-08 


16.38 


11.00 


9.178E-08 


7.770E-09 


8.47 


1.417E-08 


15.44 


13.00 


2.380E-08 


3.856E-09 


16.20 


3.816E-09 


16.03 



TABLE XVI: 



spectrum for combined centralities: 60-92% 





Pt 


Yield 


Stat. Error 


% 


Sys. Error 


% 


1.25 


9.037E-02 


4.484E-04 


0.50 


9.891E-03 


10.94 


1.75 


1.571E-02 


9.400E-05 


0.60 


1.680E-03 


10.69 


2.25 


3.397E-03 


2.836E-05 


0.83 


3.765E-04 


11.08 


2.75 


8.712E-04 


1.041E-05 


1.19 


1.025E-04 


11.76 


3.25 


2.808E-04 


4.746E-06 


1.69 


3.442E-05 


12.26 


3.75 


9.236E-05 


2.390E-06 


2.59 


1.073E-05 


11.61 


4.25 


3.461E-05 


1.282E-06 


3.70 


4.593E-06 


13.27 


4.75 


1.369E-05 


7.188E-07 


5.25 


1.857E-06 


13.57 


5.25 


6.198E-06 


4.761E-07 


7.68 


9.308E-07 


15.02 


5.75 


2.683E-06 


2.173E-07 


8.10 


4.415E-07 


16.46 


6.25 


1.514E-06 


1.486E-07 


9.81 


2.414E-07 


15.94 


6.75 


8.605E-07 


1.046E-07 


12.16 


1.589E-07 


18.46 


7.25 


4.305E-07 


6.822E-08 


15.85 


7.426E-08 


17.25 


7.75 


1.502E-07 


4.381E-08 


29.17 


2.719E-08 


18.10 


8.50 


9.326E-08 


2.035E-08 


21.82 


1.706E-08 


18.29 


9.50 


5.374E-08 


1.486E-08 


27.65 


9.997E-09 


18.60 


11.00 


1.101E-08 


4.162E-09 


37.80 


2.062E-09 


18.73 


13.00 


1.478E-09 


1.478E-09 


100.00 


2.883E-10 


19.51 
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